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ABSTRACT 

The Sloan Digital Sky Survey (SDSS) provides data on several hundred thousand galax- 
ies. Precise location of these galaxies in the sky, along with information about their lumi- 
nosities and line-of-sight (Doppler) velocities allows one to construct a three-dimensional 
map of their location and estimate their linc-of-sight velocity dispersion. This informa- 
tion, in principle, allows one to test dynamical gravity models, specifically models of 
satellite galaxy velocity dispersions near massive hosts. A key difficulty is the separation 
of true satellites from interlopers. We sidestep this problem by not attempting to derive 
satellite galaxy velocity dispersions from the data, but instead incorporate an interloper 
background into the mathematical models and compare the result to the actual data. 
We find that due to the presence of interlopers, it is not possible to exclude several 
gravitational theories on the basis of the SDSS data. 
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1 INTRODUCTION 



Recently, IKlvpin fc Pradal (|20Qgh presented an analysis of 
galaxy observations of the Sloan Digital Sky Survey (SDSS, 
http : / /www . sds s . org/] ) to test gravity and dark matter in 
the peripheral parts of galaxies at distances 50-400 kpc from 
the centers of galaxies. This field of extragalactic astronomy 
provide s one of the main arguments for the presen ce of dark 
matter (|Zaritskv fc Whitdfl994l; iPrada et ai1l2003h . 

The analysis of lKlypin fc Pradal (|2009h be gins with iden- 
tifying candidate host galaxies and candidate satellite galax- 
ies based on their relative positions in the sky, relative veloc- 
ities, and relative luminosities. After a candidate population 
of hosts and satellites has been identified, an ad-hoc math- 
ematical model is used to separate the (presumed constant) 
background of interlopers from actual satellites. This mathe- 
matical model yields a velocity dispersion profile for the pre- 
sumed satellites that is then checked against theory. 

In the present work we propose an alternative approach 
that altogether avoids the difficult issue of identifying satel- 
lites vs. interlopers. Rather than attempting to subtract the 
interloper population from the data in order to construct a 
data set that is then hoped to represent the satellite popula- 
tion correctly, we endeavor to model the actual data instead, 
by adding an interloper population to the velocity dispersions 
predicted by various gravity theories. Crudely put, we extend 
the theory to model the data correctly, rather than massaging 
the data to fit within the constraints of a limited model. 

In the first section of our paper, we offer a detailed de- 
scription of our data analysis. In the second part, we model 
the data using three gravity theories. In addition to New- 
tonian gravity and Modified Newtonian Dynamics (MOND, 



iMilgroml |l984l ). of particul ar inter e st to u s is our Modi- 
fied G ravity Theory (MOG, iMoffafj (|2006l ): iMoffat fc Tothl 
(2009)), which has been used successfully in the past to ex - 



plain galaxy rotation curves 
galaxy cluster mass profiles 



Brownstein fc M offat 
Brownstein fc MoffatJ I2OO6J ). 



cosmological observations (jMoffat fc Tothl 120071). and grav- 
itatio nal lensing in the Bullet Cluster ([Brownstein fc Moffad 
2007) without assuming the presence of nonbaryonic dark 
matter. In the third section, we combine our satellite veloc- 
ity dispersion predictions with the observed interloper back- 
ground, and contrast the resulting predic tions, as well as 
the co ld dark matter (CDM) prediction of IKlvpin fc Pradal 
(|2009h with the SDSS data. Our conclusion is that the SDSS 
galaxy data cannot be used to exclude any of these gravi- 
tational theories, not unless an independent, nonstatistical 
method is found that can be used reliably to identify individ- 
ual interlopers. 



2 DATA ANALYSIS 

The SDSS3 Data Release 6 (DR6) provides imaging data over 
9500 deg 2 in five photometric bands. Galaxy spectra are de- 
termined by CCD imaging and t he SDSS 2.5 m telescope 
on A pache Point, New Mexico (|Adelman-McCarthv et all 
120081 ). Over half a million galaxies brighter than Petrosian 
r-magnitude 17.77 over 7400 deg 2 are included in the SDSS 
data with a redshift accuracy better than 30 km/s. 

Due to the complications of calculating modified gravity 
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Figure 1. Line-of-sight velocities (top left) and number densities (bottom left) of candidate satellite galaxies (full sample) as a function of 
projected distance from the candidate host. After removal of candidate interlopers, the number density (bottom right) follows a power-law 
profile with an exponent of — 1.5 (top right). 



for non-spherical objects, iKlypin fc Pradal (|2009l ) restricted 
their analysis only to red galaxies, the vast majority of which 
are either elliptical galaxies or are dominated by bulges. We 
followed a similar strategy, restricting our selection of candi- 
date host galaxies to isolated red galaxies. We also restricted 
our selection to galaxies with a recession velocity between 
3000 km/s and 25000 km/s, which yielded approximately 
234,000 galaxies in total. 

We began our analysis by obtaining a data set from the 
SDSS. We obtained sky positions, spectra, and extinction- 
corrected magnitudes for 687,423 galaxies. We adjusted the 
data set by accounting for the motion of the solar system. We 
then processed the result using a C-language program that se- 
lected, as candidate hosts, isolated red galaxies with no other 
galaxy within a projected distance of 1 Mpc and a luminosity 
more than 25% that of the candidate host. We then identi- 
fied as candidate satellites galaxies that were within 1 Mpc of 
projected distance from the candidate host, and had a line- 
of-sight redshift velocity of less than 1,500 km/s relative to 
the candidate host. These candidate satellites were binned 
by distance. The computation yielded 3,589 hosts with 8,156 
satellites. Of these, 121 satellites (or about 1.5% of the to- 
tal) were assigned to multiple hosts; no attempt was made to 
eliminate these duplicates. 

The number density of the candidate satellites (Figure [1] 
bottom left) suggests that many of these galaxies are not, in 
fact, satellites. Indeed, if dim galaxies were distributed com- 



pletely randomly, with no relation to the candidate host, we 
would expect a number density that increases linearly with 
projected radius. The actual number density plot appears to 
be a distribution with a peak at ~100 kpc, superimposed on 
just such a linear density profile. Subtracting the linear den- 
sity profile yields the plot in Figure[T] (bottom right), which is 
a power law profile with exponent —1.5, as shown in FigureQ] 
(top right). (This corresponds to a parameter of 7 ~ —2.5 in 
the Jeans equation, discussed below.) 

3 SATELLITE GALAXY VELOCITY 
DISPERSION 

Predictions for modified gravity can be made by solving the 
Jeans equation, which gives the ra dial velocity dispersion 
af(r) as a function of radial distance. IKlypin fc Pradal (2009) 
find that neither Newtonian gravity without nonbaryonic 
dark matter, nor Modified New tonian Dynamics (M OND) are 
compatible with observations. [Angus et all (|2008l ). however, 
demonstrated that a suitably chosen anisotropic model and 
appropriately chosen galaxy masses can be used to achieve a 
good fit for MOND. 

Radial velocity dispersions in a spherically symmetric 
gravitational field can be c omputed using the Jeans equation 
(|Binnev fc Tremaindll987r i: 

d(ya 2 r ) 2v 2 d$ 

— h — /3a r = -v^-, 1) 

ar r ar 
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where v is the spatial number density of particles, v r is the 
radial velocity, f3[r) = 1— [<r 2 (r)+cx5(r)]/2cr 2 (r) is the velocity 
anisotropy, <3>(r) is the gravitational potential, and we are 
using spherical coordinates r, 9, <f>. We can write Eq.Q in 
the form 



do 2 . Aa 2 , . 

-rr + —r = -g(r ■ 

Here, we have 
A = 2/3(r)+ 7 (r), 



(2) 



(3) 



where 7(7-) = d\nv(r) / 'dlnr. 

If we assume that the velocity distribution of satellite 
galaxies is isotropic, j3 = 0. In general, f3 needs to be neither 
zero nor constant. The number density of candidate satellites 
favors a value of 7 ~ —2.5, corresponding to the observed 
power-law surface density with exponent —1.5. 

The observed velocity dispersion is along the observer's 
line of sight, seen as a function of the projected distance from 
the host galaxy. Therefore, it is necessary to integrate veloc- 
ities along the line-of-sight: 
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Figure 2. SDSS velocity dispersion data vs. gravity theory model 
predictions. The data (black dots with error bars) includes host 
galaxies in the —20.5 > M g > —21.6 luminosity range. Modified 
Gravity (MOG) is shown as the vertical shaded (blue) line. MOND 
is the horizontal shaded (red) area, mostly hidden by the MOG 
prediction. The Newtonian prediction without nonbaryonic dark 
matter is shown as a solid (green) area, while the CDM prediction is 
the dotted grey area. Similar results are obtained for the luminosity 
-21.1 and 



f[y 2 + (i 

olos(-R) = 



P)R 2 



2 a 2 (y)v(y) dy 



(4) 



where v is the spatial number density of satellite galaxies as 
a function of distance from the host galaxy, and y is related 
to the projected distance R and 3-dimensional distance r by 



r 2 = R 2 + y 2 . 



(5) 



Changing integration variables to eliminate y, we can express 
the observed line-of-sight velocity dispersion as a function of 
projected distance as 



l^ 2 



- /3R 2 )a 2 (r)u(r)/rVr 2 - R 2 dr 



(6) 



J rv(r) I \]r 2 — R? dr 



From the field equations derived from the MOG action, 
we obtain the modi fied Newtonian acceleration law for weak 
gravitational fields l|Moffatl [20061 : iMoffat fc Tothll2009l ) of a 
point source with mass M: 



J/mogW = {1 + a [1 - e"* r (1 + fir)] } . 



(7) 



where Gn is the Newtonian gravitational constant, while the 
MOG parameters a and fi determine the coupling strength 
of the "fifth force" vector <^> M to baryon matter and the range 

of the force, respect ively. 

In recent work jMoffat fc Tothll2009l ). we have been able 
to develop formulae that predict the values of the a and fi 
parameters, in the form 



D_ 



M 



(VM + E) 2 
where the parameters 
Goo — 20Gjv, 



Gn 



D 
E 



6250 Mq /2 



kpc 



25000 M, 



1/2 



(8) 
(9) 

(10) 

(11) 
(12) 



are determine d from galaxy rotatio n curves and cosmological 
observations (|Moffat fc Tothll2009l ). 

The MOND acceleration (7mond is given by the solution 
of the non-linear equation 



JMONDjK 



^ lffM< -> 



G N M(r) 



(13) 



where M is the mass of only baryons and Oo = 1.2 x 
10 _8 cm sec -2 . The fo rm of the fun ction n(x) originally 
proposed by Milgrom jMilgroml fl984h is given by fj,(x) = 
x/y/l + x 2 ; however, better fits and bette r asymptotic be- 
havio r are achieved using fi(x) — x / (1 + x) (jKlypin fc Pradal 
l2009h . 

Following in the footsteps o f lKlvpin fc Pradal j2009), we 
grouped satellite galaxy velocities for host galaxies in two 

luminosity ranges: -20.5 > Mg 1 ' > -21.1, and -21.1 > 

(2) 

Mg > —21.6. T he corresponding masses for the host galax- 



ies, calc ulated by Klypin fc Pradal |2009l) on the basis of the 
work of lBell fc de Jond (|200ll ). are 



Mi 



(i) 



Ml 



(2) 



7.2 x 10 iu M @ (-20.5 > M ( g L> > -21.1), (14) 
1.5 x 10 11 M© (-21.1 > M {2) > -21.6). (15) 



We used these values to obtain two sets of predictions 
for each theory, using ft = 0, 7 = —2.5. 



4 THE INTERLOPER BACKGROUND 

Having obtained the velocity dispersion for satellite galax- 
ies around a host galaxy, we now turn our attention to the 
interloper population. 

The actual data consists of host galaxies, satellites, and 
an effectively random interloper background. When satellites 
and interlopers are binned by projected distance from host 
galaxies, the result can be modeled symbolically as 



N(R) ± 5N(R) = N mt {R) + N int (R), 



(16) 



where N(R) is the number of galaxies in the bin at projected 
radius R, N sa t(R) and Mnt(-R) are the number of satellites 
and interlopers, respectively, in that same bin, and S is used 
to represent the sampling error. 
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This is not the approach taken by iKlypin fc Pradal 
(2009), however. Instead, they elected to subtract a modeled 
interloper background from the observed number density of 
satellites, and then compare that to a model representing only 
satellite galaxies. In effect, they used 



[N(R) - N int (R)] ± tfAU(i?) = N mt (R). 



(17) 



Assuming that the sampling error of satellites and interlopers 
are independent, we have 



SN(R) = y/SNL(R) + SN? nt (R) > SN^(R), 



(18) 



leading to potentially misleading conclusions about the extent 
to which the galaxy sample can be used to constrain alternate 
gravity models. It was this realiza tion that led us t o rep eat 
some of the analysis performed bv lKlypin fc Pra da (2009). 

For this reason, in our analysis we endeavor to model 
the actual observation, by estimating both satellite galaxy 
velocity dispersions in accordance with the previous section 
and the velocity dispersion of the interloper background. We 
assume a constant (i.e., independent of distance or sky posi- 
tion) interloper background. We find that in a data bin of size 
Arint = 25 kpc between 975 kpc and 1 Mpc) a distance at 
which no true satellites are expected to exist, the interloper 
velocity dispersion is <Tmt ~ 622.5 km/s. The average line- 
of-sight velocity of galaxies in this bin is 29.9 km/s, which 
agrees well with the presumed accuracy of SDSS line-of-sight 
velocity estimates; therefore, we use 29.9 km/s as the stan- 
dard deviation of the velocity dispersion data. We note the 
number of interlopers found in this bin, ./Vint = 152. Consid- 
ering now a predicted satellite galaxy velocity dispersion in 
a bin of size Ar at projected radius r, the combined velocity 
dispersion can be calculated as 

2 2 I "-, >"sat A r 8a t iVint \ , 2 r sat 



(19) 



V n»t Annt N J r mt Ar int N 

where N is the observed number of candidate satellites in the 
selected bin. 

In Figure [2] we plotted theoretical predictions — that is, 
solutions of Eq. ((2j using the Newtonian, MOG JJ), and 
MOND (|13|) accelerations for the host galaxy mass range de- 
fined by Eqs. (|14[) and (| 15p — against the binned SDSS data 
from Figure [T] (lower left). The data points are taken from 
the SDSS data set, with host galaxy luminosities M g limited 
to —20.5 > M g > —21.6. Four theoretical predictions are also 
shown for this luminosity range. The predictions for Newto- 
nian theory without nonbaryonic dark matter, MOND, and 
MOG were calculated using the approach described in this 
pape r, while the standard m odel CDM prediction is taken 
from lKlypin fc Pradal (|2009h . 

As the figure shows, the theories are practically indis- 
tinguishable (indeed, they are barely distinguishable from a 
"null theory" of stationary satellite galaxies!) This is a result 
of the fact that in the observed velocity dispersion, interlop- 
ers dominate due to their much larger velocities, even when 
satellites dominate in numbers. As an example, at 125 kpc, 
for hosts in the —21.1 > M g > —21.6 luminosity bin, the 
CDM prediction is 219.8 km/s, while the Newtonian pre- 
diction without nonbaryonic dark matter is only 31.8 km/s. 
However, once the interloper background is considered with 
a velocity dispersion of 622.5 km/s, using an estimated satel- 
lite number density of 140.7 and interloper number density 
of 17.3 in this data bin, the resulting velocity dispersions are 
292.3 km/s for CDM and 208.2 km/s for the Newtonian the- 
ory without nonbaryonic dark matter; a much less dramatic 
difference. 



For this reason, it seems futile to use this type of statisti- 
cal analysis of satellite galaxies to distinguish between CDM 
models on the one hand, and various gravitational theories on 
the other. The data appear to be consistent with CDM and 
also with MOND, MOG, and marginally, even with the New- 
tonian theory without assuming the presence of nonbaryonic 
dark matter. 



5 CONCLUSIONS 



Obser vational data p resented by Adclman-Mc Carthv et al.l 
(120081 ) and studied bv lKlypin fc Pradal l|2009h are viewed as 
evidence of the success of the ACDM model. The usual ap- 
proach relies on the critical step of interloper removal, before 
the data is compared against predictions. We argue that this 
approach is fundamentally flawed: rather than attempting to 
remove interlopers from the data, we must add the interloper 
background theoretical predictions, in order to predict obser- 
vational values. When we carry out this approach, we find the 
CDM and gravity theory predictions cannot be distinguished 
and that the SDSS data cannot be used to falsify any of the 
theories we examined. 



ACKNOWLEDGMENTS 

The research was partially supported by National Research 
Council of Canada. Research at the Perimeter Institute 
for Theoretical Physics is supported by the Government of 
Canada through NSERC and by the Province of Ontario 
through the Ministry of Research and Innovation (MRI). 



REFERENCES 

Adelman-McCarthy, J. K., Agueros, M. A., Allam, S. S., 
Allende Prieto, C, Anderson, K. S. J., Anderson, S. F., 
Annis, J., Bahcall, N. A., Bailer- Jones, C. A. L., Baldry, 
I. K., et al. 2008, ApJS, 175, 297 

Angus, G. W., Famaey, B., Tiret, O., Combes, F., & Zhao, 
H. S. 2008, Mon. Not. R. Astron. Soc, 383, LI 

Bell, E. F., & de Jong, R. S. 2001, Astrophys. J., 550, 212 

Binney, J., & Tremaine, S. 1987, Galactic Dynamics (Prince- 
ton University Press) 

Brownstein, J. R., & Moffat, J. W. 2006a, Mon. Not. R. 
Astron. Soc, 367, 527 

— . 2006b, Astrophys. J., 636, 721 

— . 2007, Mon. Not. R. Astron. Soc, 382 (1), 29 

Klypin, A., & Prada, F. 2009, Astrophys. J., 690, 1488 

Milgrom, M. 1984, Astrophys. J., 270, 365 

Moffat, J. W. 2006, Journal of Cosmology and Astroparticle 
Physics, 2006, 004 

Moffat, J. W., & Toth, V. T. 2007, ArXiv, 0710.0364 [astro- 
ph] 

— . 2009, Class. Quant. Grav., 26, 085002 

Prada, F., Vitvitska, M., Klypin, A., Holtzman, J. A., 
Schlegel, D. J., Grebel, E. K., Rix, H.-W., Brinkmann, J., 
McKay, T. A., & Csabai, I. 2003, Astrophys. J., 598, 260 

Zaritsky, D., & White, S. D. M. 1994, Astrophys. J., 435, 
599 



© 0000 RAS, MNRAS 000, 000-000 



